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Quantum defect theory is applied to (time-dependent) density-functional calculations of Rydberg 
series for closed shell atoms: He, Be, and Ne. The performance and behavior of such calculations 
is much better quantified and understood in terms of the quantum defect, rather than transition 
energies. 



I. INTRODUCTION 

Time-dependent density functional theory 
(TDDFT) [l| has enjoyed a recent surge in popu- 
larity for calculating excited-state energies of atoms, 
molecules, clusters, and solids 0, 0- TDDFT has fea- 
tures similar to ground-state density functional theory 
(DFT) : It produces useful accuracy at a fraction of the 
computational cost of ab initio methods Q , but reliabil- 
ity depends on the approximate functionals used As 
more and more practictioners in many subfields of com- 
putational science use TDDFT, there are ever increasing 
numbers of implementations. Calculations on atoms and 
sequences of atoms are often used to benchmark new 
implementations, or to show how well TDDFT works 
in the simplest cases. These are particularly useful, as 
much highly accurate data, both from experiment and 
accurate wave function calculations, are available for 
these systems. 

A well-known difficulty that hampered even the earli- 
est calculations of excitations in atoms with TDDFT 5] 
is the incorrect asymptotic behavior of the ground-state 
potentials of common density functional approximations. 
It has long been known that the exact KS potential for 
a closed shell atom decays as — 1/r at large r, where 
r is the distance from the nucleus 

H 0- Typical ap- 
proximations for the ground-state, such as the local den- 
sity approximation (LDA) , gradient corrected functionals 
(GGAs), and hybrid functionals (see for example Ref. 8] 
and references therein), have potentials that decay too 
rapidly with r. Only those with the correct asymptotic 
behavior support the Rydberg series of transitions, an in- 
finite number of transitions that merge with the contin- 
uum at the ionization threshold. Orbital-dependent func- 
tionals capture this behavior naturally, as does the Van 
Leeuwen-Baerends potential approximation (LB94) 0, 
which was designed to have an asymptotically correct 
behavior. Several recent methods have been suggested 
for correcting the standard functionals to produce the 
long-ranged tail [lfj. Illl fl2l|. In any event, our present 
work applies only to long-ranged potentials. 

We argue here that the way in which results have been 
calculated and reported for atoms is far from optimal. We 
show that long lists of transition frequencies for Rydberg 
series converging to the ionization threshold are not the 
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FIG. 1: Energy level diagram for the helium atom. 



best way to report such calculations. Instead, the well- 
developed theory of the quantum defect [H|, used for 
decades in atomic physics, is ideal for this purpose. We 
show that, for each Rydberg series, i.e., for each value of 
angular momentum two or three numbers completely 
characterize all the information in the infinite series. Fur- 
thermore, the quantum defect is a much more demand- 
ing test of excitation energies, and methods that appear 
to have only small energetic errors can yield quite poor 
quantum defect behavior. Also, shifts in orbital energies, 
such as the missing correlation contribution in the ex- 
act exchange HOMO energy, have no effect on the quan- 
tum defect, so that the quality of a potential can be as- 
sessed without being influenced by such errors. We also 
find , that the quantum defect of the exact ground- 
state Kohn-Sham (KS) potential, is sometimes, but not 
always, a good starting point for approximations to the 
true quantum defect. 
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FIG. 2: The exact s KS quantum defect and the exact singlet 
and triplet quantum defects fl^l of He and their parabolic fits. 



II. THEORY 

In Fig. ^ we show the orbital energy level diagram of 
the helium atom. The zero is set at the onset of the 
continuum, marked with a dotted line. For closed shell 
atoms and for any spherical one-electron potential that 
decays as —1/r at large distances, the bound-state tran- 
sitions form a Rydberg series with frequencies: 



1 



2(n - n n i) 2 



(1) 



where I is the ionization potential, and [i n i is called the 
quantum defect. We use atomic units (e 2 = H = m e = 1) 
throughout. The value of the quantum defect is that, for 
real atoms, quantum defects depend only weakly on the 
principle quantum number n for large n and converge 
to a finite value in the limit n — ^ oo. In fact, accord- 
ing to Seaton's theorem 0], the quantum defect is a 
smooth function of energy as E — > 0, and merges con- 
tinuously with the phase-shift (relative to pure Coulomb 
scattering) divided by tt. In TableQ] we report extremely 
accurate results from wavefunction calculations for the 
helium atom. Wc show singlet and triplet values that 
have been obtained by Drake [l6|. We also give results 
from the exact ground-state KS potential, as found by 
Umrigar and Gonze |17| . We say more on how we ob- 
tained the KS values in section ITlI Al On the left are the 
transition frequencies, while on the right are the corre- 
sponding quantum defects. Note how small the differ- 
ences between transitions become as one climbs up the 
ladder, and yet the quantum defect remains finite and 
converges to a definite value. 

In Fig. |21 we show the exact s KS quantum defect and 
the singlet and triplet quantum defects corresponding to 
accurate wave function results for helium as sym- 
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FIG. 3: Same as Figure [5] plotted against energy. 

bols. This is the way the defect was plotted in Ref. 
for the case of the Ne atom. However, it is more appro- 
priate still to plot the defect as a function of energy, as 
shown in Fig. |3] This clearly illustrates that the quan- 
tum defect is a smooth function of energy, and will be 
well approximated a polynomial of some low order p: 



i=0 



E = uj-I. 



(2) 



We choose to optimize the fit over the entire range of ex- 
itation energies, not just about E = (the [ii are simply 
related to the a, b, and c coefficients in Ref. For 
example the data in Fig. |3| can be accurately described 
by a straight line. We report these coefficients in Ta- 
ble [H] and give the quantum defects obtained from these 
coefficients as continuous lines in Figs. |2] and [3J We also 
added the transition energies for He obtained from the 
KS fit to Tabled The maximun error in the fit is 3/100 
mH so it is essentially exact for all purposes of this pa- 
per. We shall see that in other cases just two or even 
three coefficients are not enough to describe the data ac- 
curately. In those cases we will need to decide when to 
stop adding more coefficients to fit the data, since we 
would like to describe the data with as few coefficients 
as possible. Therefore we will look at the maximum ab- 
solute error in the quantum defect. By which we mean 
that we recalculate the quantum defects from the coef- 
ficients and look at the absolute differences between the 
fitted and original values, and then look at the largest dif- 
ference. We will stop adding more coefficients once the 
maximum absolute error is smaller than 0.001, or, in case 
this value is not reached, whenever the absolute errors do 
not change much when adding more coefficients. We give 
the value for the error also in Table [H] Any approximate 
ground-state KS potential suggested for use in TDDFT 
should have its coefficients compared with the KS num- 
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TABLE I: 


Transition energies 


for He atom [a.u.]. 


The ionization 


energies 


are 0.9037244 a.u. 


for the singlet and triplet case 


ana u.yuo / z 


a.u. for the KS potential. 
















AE 








Quantum defect 




Transition 


Singlet a 


Triplet ° 


KS 6 


KS fit 


Singlet 


Triplet 


KS 


Is — > 2s 


0.7577503 


0.7284949 


0.74599 


0.74596 


0.1492525 


0.3107982 


0.21957 


Is — > 3s 


0.8424524 


0.8350353 


0.83917 


0.83920 


0.1433699 


0.3020047 


0.21689 


Is -» 4s 


0.8701376 


0.8672122 


0.86883 


0.86882 


0.1416535 


0.2994464 


0.21492 


Is — > 5s 


0.8825475 


0.8883469 


0.88189 


0.88189 


0.1409171 


0.2983583 


0.21459 


Is — > 6s 


0.8891613 


0.8925944 


0.88879 


0.88879 


0.1405330 


0.2977954 


0.21441 


Is -» 7s 


0.8930986 


0.8925945 


0.89287 


0.89287 


0.1403072 


0.2974666 


0.21429 



"Accurate non-relativistic calculations from Ref. I Hi . 
'The differences betwee n th e KS eigenvalues obtained with the 
exact potential from Ref. Il7l . 



TABLE II: Fit coefficients fio, and [ii for the s quantum 
defect of He. By 'Max. AE' we mean the maximum absolute 
error as explained in the text. 





Singlet" 


KS b 


Triplet" 


Ho 


0.1396 


0.2139 


0.2965 


Mi 


-0.0655 


-0.0370 


-0.0811 


Max. AE 


0.0002 


0.0006 


0.0001 



"Obtained from non-relativistic calculations of the orbital energies 
from Ref. 1773 . 

'KS values obtained with the exact potential from Ref. |l7j . 



bers in the table, while any approximate xc-kernel should 
have its coefficients compared with the singlet and triplet 
case. 



TABLE III: The fii from the exact ground-state KS potential 
for He, Be, and Ne. 





He a 


Be" 


Ne 6 


s 


I 


0.9037 


0.3426 


0.7945 


/to 


0.2139 


0.7164 


1.3125 


Ml 


-0.0370 


-0.2185 


-0.1807 


Max. AE 


0.0006 


0.0002 


0.0007 


P 




0.0164 


0.3587 


0.8304 


Mi 


0.0289 


-0.3377 


-0.3500 






0.6112 




Max. AE 


0.0009 


0.0003 


0.001 



Ref. IT' 



III. GROUND-STATE KOHN-SHAM 
POTENTIALS 

All linear response TDDFT calculations of excitations 
begin from the occupied to unoccupied transitions of the 
ground-state KS potential. In this section, we analyse 
some of the most popular approximations using quantum 
defect theory. 



A. Computational details 

All ground-state DFT results shown here are calcu- 
lated with a modified OEP (optimized potential model) 
program 0, 0, |2(j. This program is basis set inde- 
pendent, works with a radial grid, and both the ener- 
gies and the potentials are optimized in a self-consistent 
way. The exact-exchange (x-only) OEP is already in- 
cluded in this program. We also did calculations with the 
LB94 ;9j potential by Van Leeuwen and Baerends, which 
was not available in the program. We implemented this 
functional by adding the LB94 correction to the LDA 
xc-potential and let the program optimize this potential 
in a self consistent manner. The program has the abil- 
ity to read in the accura te p otentials by Umrigar et al. 
(He [13, Be |2ll, and Ne [2lf). These accurate potentials 



were only known by us up to a particular radius. In or- 
der to calculate energies for higher n values we needed to 
increase this radius. We did this by adding a — 1/r tail 
to u xc and a Z/r tail to v rm H and we checked that the 
transition was smooth. We also made sure our values are 
converged with the number of gridpoints. When we used 
the accurate potentials, we did not allow the program to 
self consistently change the potential. 

The maximum n value for which we could still do very 
accurate calculations is n = 7 for He and n = 9 for Be 
and Ne. 



B. Exact results 

As we have seen the numbers, /i^, contain all the in- 
formation needed to characterize a given Rydberg series, 
and make tables of the actual transition frequencies re- 
dundant. In Table IIIII we report the coefficients for s 
and p KS quantum defects for the He, Be and Ne atoms 
obtained with accurate xc potentials. We use the transi- 
tions up to n — 9 for the fit in case of Be and Ne, and up 
to n = 7 in case of He. 

When we compare the s (I = 0) values with p (I = 1), 
the asymptotic KS quantum defect is smaller for p in all 
cases. This reflects the lesser importance of the inner part 
of the KS potential relative to the angular momentum 
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FIG. 4: The exact, LB94, and x-only potentials for beryllium 
together with — 1/r. 



barrier as I grows. However, the curvature of the s and 
p quantum defect is similar. 



C. Approximations 

In this section, we demonstrate our methodology by 
testing two common approximations for the ground-state 
KS potential. These are exact exchange OEP [20j and 
LB94 9j ■ Exact exchange calculations are more demand- 
ing than traditional DFT calculations, but are becom- 
ing pop ular because of the high quality of the poten- 
tial \'2\j . |23| . On the other hand, LB94 provides an 
asymptotically correct potential at little extra cost be- 
yond traditional DFT [H H H|. Fi S- H shows both 
these potentials for the Be atom together with the exact 
potential [2l| . Our figure for the exact potential is differ- 
ent from the one in Ref. 0, since in that reference they 
use a different accurate density to derive the potential. 
We give our values for the fit parameters in Table llVl In 
Fig. we show the p Be quantum defect obtained with 
LB94, OEP, and KS and we show the fit as continuous 
lines. Fig. [S] immediately shows the high quality of the 
OEP potential. The quantum defect curve is almost iden- 
tical to the exact one, being offset by about 0.1 (see Ta- 
ble [Wj . On the other hand the quantum defect of LB94 
is poor, and this is true for all cases studied. From Fig.0] 
we can see that the OEP is much closer to the true po- 
tential than LB94 and also approaches — 1/r faster. This 
shows that just having a potential that is asymptotically 
correct is not enough to get a good quantum defect. 

Another thing we have not mentioned so far is that 
a potential that gives a wrong ionization potential does 
not necessarily give a bad quantum defect. The origin of 
this is that the quantum defect is obtained from lu„i — I 
and can take a shift of energy levels into account. Con- 
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FIG. 5: The Be p quantum defect of LB94, OEP, and KS, 
and their best fits. 



sider the exact exchange results for the atoms. Typically, 
these quantum defects are accurate to 0.1. Thus, using 
the exact ionization potential with the exchange quan- 
tum defects, yields highly accurate transition frequen- 
cies, i.e., the most significant error in OEP excitations 
is due to the missing correlation contribution to the po- 
sition of the HOMO. On the other hand, we see that 
LB94, while asymptotically correct and sometimes hav- 
ing a highly accurate ionization potential, has much less 
accurate quantum defects. 



IV. 



TDDFT RESULTS 



In the previous sections we saw that the KS quantum 
defects are typically lying in between the exact singlet 
and triplet quantum defects In order to calculate 

these singlet and triplet quantum defects within DFT 
the usual method of choice is TDDFT (within the linear 
response regime). Apart from a ground-state potential 
one needs to choose an xc-kernel as well. In this section 
we obtain quantum defects from excitation energies ob- 
tained with different xc-kernels and ground-state poten- 
tials. We focus mainly on the Be atom, but we also give 
expansion coefficients for He and Ne quantum defects. 



A. Performance of TDDFT 

We concentrate on the quantum defect obtained from 
the ALDA kernel applied to the exact ground-state KS 
potential. In Table El we show s — > s and s — > p exci- 
tation energies for the Be atom. The ALDA results in 
column three are obtained by Van Gisbergen et. al. . 
For their ground-state calculations, they used the accu- 
rate potential by Umrigar and Gonze |2ll | and for the 
xc-kernel they used the ALDA. The calculations were 
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TABLE IV: The fj,i from different ground-state potentials for He, Be, and Ne. The ionization energies are in [a.u.] and are 
not included in the maximum error. 









s 






P 








OEP 


LB94 


exact 


OEP 


LB94 


exact 


He 


I 


0.9182 


0.8513 


0.9037 










Mo 


0.2170 


-0.3607 


0.2139 


0.0238 


-0.6015 


0.0164 




/'i 


-0.0401 


-2.7408 


-0.0370 


0.0319 


-2.2013 


0.0289 








-3.1852 






-1.5723 






Max. AE 


0.0000 


0.003 


0.0006 


0.0000 


0.002 


0.0009 


Be 


I 


0.3093 


0.3205 


0.3426 










Mo 


0.6623 


0.3497 


0.7164 


0.2786 


0.0115 


0.3287 




Mi 


-0.1557 


-0.7203 


-0.2185 


-0.1070 


0.4836 


-0.3377 




H2 




31.5313 




0.8453 


66.0517 


0.6112 




M3 




73.1572 






268.6094 






Max. AE 


0.0002 


0.001 


0.0002 


0.0002 


0.002 


0.0003 


Ne 


/ 


0.8508 


0.7821 


0.7945 










Mo 


1.3433 


0.8174 


1.3125 


0.8631 


0.3196 


0.8340 




Ml 


-0.2285 


-3.2880 


-0.1807 


-0.3972 


-2.8338 


-0.3500 








-11.0968 






-1.3271 






M3 




-39.2359 












Max. AE 


0.0004 


0.004 


0.0007 


0.0008 


0.003 


0.001 



TABLE V: Transition energies for the beryllium atom. 



Transition 


Expt. a 


KS b 


ALDA C 


Truncated- 
ALDA d 


Truncated- 
Hybrid 6 


WY f 


AC-LDA 9 


2s -> 3s 


0.24913 


0.24437 


0.2495 


0.2515 


0.2510 


0.2512 


0.2239 


2s 4s 


0.29728 


0.29586 


0.2977 


0.2984 


0.2985 


0.3033 


0.2680 


2s -> 5s 


0.31586 


0.31526 


0.3160 


0.3164 


0.3165 






2s -> 6s 


0.32496 


0.32466 




0.3252 


0.3254 






2s -> 2p 


0.19394 


0.13273 


0.1868 


0.1889 


0.1427 


0.1861 


0.1781 


2s -> 3p 


0.27423 


0.26937 


0.2710 


0.2714 


0.2736 


0.2749 


0.2457 


2s -> 4p 


0.30543 


0.30461 


0.3048 


0.3049 


0.3059 


0.3107 


0.2754 


2s -> 5p 


0.31949 


0.31931 




0.3194 


0.3199 






2s -> 6p 


0.32690 


0.32686 




0.3269 


0.3272 






s MAE [mH] 






0.312 


1.069 


0.943 


4.046 


27.25 


p MAE [mH] 






3.669 


1.700 


10.61 


4.592 


24.81 


Total MAE [mH] 






1.991 


1.420 


6.313 


4.374 


25.79 



"Experimental values from NIST |3 

'The differences betwee n th e KS eigenvalues obtained with the 
exact potential from Ref. 12 ll 

C ALDA calculation including all bound and unbound states from 
Ref. 

d ALDA calculation including 34 unbound states from Ref. |2jj 
e Hybrid calculation including 34 unbound states for He and 38 
unbound states for Be from Ref. I'-Hj 
^ALDA calculation with WY ground-state potential from Ref. l3ll 
'Asymptotically corrected ALDA results from Ref. l30l 



done close to the basis set limit and with high numeri- 
cal integration accuracy. As can be seen from the table 
the excitation energies for the s — > s transitions are very 
close to the experimental values with a mean average 
error (MAE) of only 0.3 mH. The s — > p excitation en- 
ergies are a bit less accurate with an MAE of 3.7 mH, 
but this is still an accurate result. The fit coefficients for 
the quantum defects of these calculations are reported in 



Table[VT| For the fit of the TDDFT results we took a less 
strict constraint to determine when to stop adding more 
coefficients. We took 0.01 instead of 0.001, reflecting the 
greater error in this data. 

In Fig. we show the s and p quantum defects corre- 
sponding to these values and we compare them with the 
bare KS and experimental results. For the s quantum 
defect, the experimental curve is essentially a straight 
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TABLE VI: The fa from different singlet TDDFT values for Be. 





Ref." 


KS 6 


ALDA C 


Truncated- ALDA d 


Truncated- Hybrid 13 


WY/ALDA / 


AC-LDA 9 


s 


I 


0.3426 


0.3426 


0.3426 


0.3426 


0.3426 


0.3493 


0.3111 


fio 


0.6752 


0.7164 


0.6536 


0.6182 


0.6010 


0.6684 


0.5997 


Ml 


-0.0133 


-0.2185 


-0.2986 


-0.4294 


-0.6947 


-0.7540 




fJ.2 


1.2345 














Max. AE 


0.0003 


0.0002 


0.004 


0.002 


0.004 


0.0000 


0.006 


P 


Mo 


0.3745 


0.3287 


0.3327 


0.3427 


0.3087 


0.3484 


0.2413 


Hi 


1.1612 


-0.3377 


-1.3172 


-0.9634 


0.2505 


-1.9702 


-1.1215 


fl2 


1.5678 


0.6112 


-13.5673 


-12.4534 


3.9961 


-15.7792 


-18.6196 


H3 


21.4587 














Max. AE 


0.0002 


0.0003 


0.0000 


0.003 


0.004 


0.0000 


0.0000 



"Experimental values from NIST |2(1 

'The differences betwee n th e KS eigenvalues obtained with the 
exact potential from Ref. 12 ll 

C ALDA calculation including all bound and unbound states from 
Ref. 

d ALDA calculation including 34 unbound states from Ref. |2i3 
e Hybrid calculation including 38 unbound states from Ref. |2Sl 
■'ALDA calculation with WY ground-state potential from Ref. 30] 
'Asymptotically corrected ALDA results from Ref. l30l 
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FIG. 6: The s and p experimental, KS, and ALDA singlet 
quantum defects for Be. The continuous lines are the fi^ (E) 
fits. 



line, with a small negative slope. The p quantum defect, 
on the other hand, is much more curved, with a large 
positive slope, due to the much lower 2s — > 2p transi- 
tion. The exact KS potential has been touted as a good 
approximation to the experimental results |3lj |. This is 
clearly true for the s curves but not so for the p quantum 
defect. In that case the KS quantum defect, while being 
very close to the experimental value as n — > oo has the 
wrong behavior as a function of E. 

In both the s and p case it is clear that doing a full 



TDDFT calculation considerably improves upon the bare 
KS results. The ALDA does slightly overcorrect the s 
quantum defect, underestimating the value at E = 0. In 
case of the p quantum defect the ALDA tends to correct 
for the opposite slope of the KS values compared to the 
experimental values, but the correction is not complete, 
leading to a curved line. We also see that the s quantum 
defects as obtained with the ALDA are much better than 
the p quantum defects, even though the MAE is only a 
few mH in the last case. 



B. Truncating Casida's equation 

Another set of ALDA calculations were performed by 
Petersilka, Gross, and Burke |28j. The difference between 
their calculation and that of Van Gisbergen is that they 
truncate the summation over states in the response func- 
tion, including only poles of bound states and neglecting 
continuum contributions. They included the lowest 34 
bound states of s and p symmetry in their calculations. 
Just like in case of the full ALDA calculations the ground- 
state was determined with the potential of Umrigar and 
Gonze ,21] . We show these results in column four of Ta- 
ble[Vl Again the results are close to the experimental val- 
ues with an MAE of 1.1 mH for the s — > s transitions and 
an MAE of 1.7 mH for the s — > p transitions. In Fig. [7| 
we show the s and p quantum defects corresponding to 
these values and we compare them with the full ALDA 
results. For the s quantum defect we see that the quan- 
tum defect is slightly below the ALDA values, and the 
slope is too large, leading to a smaller asymptotic quan- 
tum defect. For the p quantum defect we see that there 
is not much difference between the truncated and the full 
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FIG. 7: The s and p ALDA and truncated- ALDA singlet FIG. 8: The s and p experimental, ALDA, and WY quantum 
quantum defects for Be. The continuous lines are the fi^ (E) defects for Be. The continuous lines are the ^ P \E) fits, 
fits. 



results. So the effect of the truncation in case of Be is 
not so great. For He the difference is relatively larger, 
because the TDDFT corrections are so small. This can 
be seen from the values in Tables fVTI and IVTT1 which we 
will discuss in section flVFI 



C. Quality of the ground-state potential 

For all TDDFT methods we have described so far, the 
ground-state was calculated with the exact potential of 
Umrigar and Gonze |2lj . This eliminates errors in the 
ground-state so one can compare the effects of using a 
different xc-kernel. But for practical calculations such 
an accurate potential is not available. This motivates the 
development of other accurate ground-state potentials to 
calculate the Rydberg series. 

Wu and Yang 12] obtained an accurate potential by 
a direct optimization method that allows them to calcu- 
late the potential from a given electronic density. This 
density is obtained from a coupled-cluster singles and 
doubles (CCSD) calculation. In Table W\ we show the ex- 
citation energies corresponding to the WY potentials and 
the error is just a few mH. In Fig. [H] we show the quan- 
tum defects obtained with this potentials and compare 
them with the ALDA and experimental values. Only the 
asymptotic value of the quantum defect is accurate. Er- 
rors due to (very small) errors in the ground-state density 
are visible for all other energies, and can be comparable 
to the xc-kernel itself. The WY method is therefore still 
very promising, but clearly requires a very accurate in- 
put density. Quantum defect analysis should prove very 
useful for testing WY-type calculations, for example, for 
comparing basis set errors with errors due to the level of 



calculation used for obtaining the input density. 



D. Testing approximate kernels 

In an attempt to improve the xc-kernel, Burke, Peter- 
silka, and Gross suggested the following form, 

f TT = f TT tU = fTXALDA /on 

J xc Jx S JXC jxc ■ \ u / 

This form is based on the fact that the parallel-spin con- 
tribution is well described in the exact exchange case, 
because of a cancellation of exchange contributions. But 
an exact exchange treatment misses the significant anti- 
parallel correlation contribution, leading to too large sin- 
glet/triplet splittings. Therefore it is recommended that 
for the anti-parallel kernel one uses the ALDA. We show 
the excitation energies obtained with this kernel in Ta- 
ble [V] It should be noted that just as in case of the 
truncated- ALDA calculations, the number of states in- 
cluded in the hybrid calculation is also limited. Namely, 
34 states in case of He and 38 states in case of Be. There- 
fore we shall denote the method by truncated-hybrid. 
The ground-state of these calculations was again calcu- 
lated with the exact potential of Umrigar and Gonze [2l| . 
Apart from the 2s — > 2p transition the truncated-hybrid 
results are very similar to the truncated- ALDA results. 
It is the large error in the 2s — > 2p transition that leads 
to the large MAE. 

In Fig. [5] we show the s and p quantum defects cor- 
responding to the truncated-hybrid results and we com- 
pare them with the truncated-ALDA and experimental 
results. From the s quantum defect plot it can again be 
seen that the truncated-hybrid results are close to the 
truncated-ALDA values, so the kernel does not improve 
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FIG. 9: The s and p experimental, truncated- ALDA, and 
truncated-hybrid singlet quantum defects for Be. The contin- 
uous lines are the ^^{E) fits. 
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FIG. 10: The s and p experimental, ALDA, and AC- ALDA 
quantum defects for Be. The continuous lines are the fi^ (E) 
fits. 



the results in this case. For the p quantum defect, the re- 
sults actually get worse with the truncated-hybrid giving 
a slope opposite the the experimental curve, so it shifts 
the KS values in the right direction but does not correct 
for the wrong slope. The hybrid kernel does not improve 
much u pon the truncated- ALDA for Be as was also found 
in Ref. 29]. It is only good for two-electron systems, for 
which it was derived. 



E. Asymptotically corrected ground-state 
potentials 

Wu, Ayers, and Yang obtained an asymptoti- 
cally corrected LDA (AC-LDA) potential by a varia- 
tional method that forces the potential to have the cor- 
rect asymptotic behavior. This is a pure DFT treatment 
that can be applied to larger molecules. We see from 
Table [V] that the AC-LDA gives a large MAE. 

In Fig. llUl we show the quantum defects obtained with 
the AC-LDA potential and compare them with the KS 
and experimental values. In both cases the AC-LDA val- 
ues are much lower than the ALDA, which was evaluated 
on the exact ground-state KS potential. The AC-LDA 
strongly underestimates the experimental s and p quan- 
tum defects. The shape of the p curve is similar to the 
ALDA curve. From this figure, it is clear that there are 
significant errors in the underlying KS potential. 



F. Coefficients for He 

In this section we give the values of the best fit for 
He as defined in Eq. |2 and as described thereafter. The 



original excitation energies from which the quantum de- 
fects and corresponding coefficients are calculated are ob- 
tained from the same sources as the Be data described 
above. In Table IVTT1 we give results for He. The sizes of 
the quantum defects are much smaller than in Be, and 
the fractional change between KS and experiment is con- 
comitantly larger. In the case of the p quantum defects, 
they can have large opposite signs. Application of ALDA 
to the exact ground-state KS potential again works well. 

In the case of He, the truncated- hybrid result improves 
upon the truncated- ALDA. As we mentioned before this 
is because the hybrid was developed for 2-electron sys- 
tems. AC-LDA and WY coefficients of He really stand 
out as behaving much differently from the other cases 
and giving large errors that cannot be removed by adding 
more coefficients. The AC-LDA and WY coefficients of 
Be did not have this problem. 



G. Triplets 

In practice, spin decomposed TDDFT is commonly 
used, and allows prediction of singlet^triplet transitions. 
In this section we discuss the coefficients of the quantum 
defect expansion obtained from triplet excitation ener- 
gies. In Table I Villi we show the fit coefficients for He 
and Be. For the WY and AC-LDA methods there is no 
triplet data available. For the triplet case, the /Lio's are 
very well reproduced in all cases, also the other coeffi- 
cients are often close. Overall the data can in most cases 
be reproduced by only two or three coefficients. 
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TABLE VII: The m from different singlet TDDFT values for He. 





Ref." 


KS 6 


ALDA C 


Truncated- ALDA d 


Truncated- Hybrid 13 


WY/ALDA / 


AC-LDA 9 


s 


I 


0.9037 


0.9037 


0.9037 


0.9037 


0.9037 


0.9026 


0.7817 


(MO 


0.1395 


0.2139 


0.1099 


0.0223 


0.0599 


0.0226 


-0.1027 


Ml 


-0.0655 


-0.0370 


-0.1388 


-0.4537 


-0.5390 


-2.4879 


-4.3920 












-12.6664 


-25.7710 




Max. AE 


0.0002 


0.0006 


0.001 


0.006 


0.007 


0.03 


0.06 


P 


Mo 


-0.0122 


0.0164 


-0.0002 


-0.0168 


-0.0355 


-0.1039 


-0.1579 


Hi 


-0.0227 


0.0289 


-0.2077 


-0.2721 


-0.2201 


-6.1973 


-2.0656 


fj> 










-107.6068 


-13.2488 




H3 










-521.8636 






Max. AE 


0.0001 


0.0009 


0.005 


0.002 


0.008 


0.007 


0.02 



"Non-relativistic calculations from Ref. lUi 

'The differences between the KS eigenvalues obtained with the 
exact potential of Ref. |17| . 

C ALDA calculation including all bound and unbound states from 
Ref. 

d ALDA calculation including 34 unbound states from Ref. |28j 
e Hybrid calculation including 34 unbound states from Ref. |29| 
■'ALDA calculation with WY ground-state potential from Ref. 30] 
'Asymptotically corrected ALDA results from Ref. l30l 



TABLE VIII: The m from different triplet TDDFT values 
for He and Be. 











Truncated- 


Truncated- 






Ref." 


ALDA 4 


ALDA C 


Hybrid'' 




He 


Po 


0.2965 


0.2719 


0.2570 


0.3171 




Ml 


-0.0811 


-0.0892 


-0.1281 






Max. AE 


0.0001 


0.004 


0.003 


0.007 


Be 


Mo 


0.7742 


0.7706 


0.7641 


0.7683 




Pi 


0.8870 


-0.3912 


-0.4015 


-0.2607 




Max. AE 


0.001 


0.0003 


0.002 


0.002 


P 


He 


Po 


0.0684 


0.0734 


0.0657 


0.05332 




Pi 


0.0457 






0.1103 




Max. AE 


0.0002 


0.003 


0.005 


0.002 


Be 


Po 


0.3620 


0.4152 


0.4101 


0.37852 




AH 


-0.4492 


0.3553 


-0.7166 


-0.6404 




P2 


1.5829 






1.5225 




Max. AE 


0.0003 


0.004 


0.006 


0.007 



"He non-relativistic calculatio ns f rom Ref. and Be and Ne 
experimental values from NIST 12(1 

b A LDA calculation including all bound and unbound states from 
Ref. |27| 

C ALDA calculation including 34 unbound states from Ref. 28] 
''Hybrid calculation including 34 unbound states for He and 38 
unbound states for Be from Ref. 1291 



describe the Rydberg series. We calculated the quantum 
defects for He, Be, and Ne with accurate xc-potentials 
and showed that the KS quantum defect lies between the 
interacting singlet and triplet values. 

We also studied approximate ground-state KS poten- 
tials, namely, the LB94 and OEP potentials. We saw 
that the OEP results are very close to the exact KS val- 
ues. The LB94 underestimates the quantum defect in all 
cases. Again the data can be fully described by a few 
coefficients. 

In the final part of the paper we calculated the quan- 
tum defects from available TDDFT excitation energies 
from the literature. We see that the quantum defect re- 
ally amplifies the error in these cases. The TDDFT data 
can be described by only a few coefficients. 

Overall we see that while having an asymptotically cor- 
rect potential guarantees the existence of a Rydberg se- 
ries, it is not necessarily a good one. This is especially 
the case for AC-LDA TDDFT and LB94 ground state 
values, and to lesser extent even for the WY TDDFT 
values. 



V. CONCLUSIONS 

In the first part of this paper we have shown that the 
quantum defect is a valuable quantity when reporting Ry- 
dberg states. A plot of the quantum defect can be more 
insightful than a list of excitation energies. The quantum 
defect can also be fitted to an expansion around E = 
and a finite number of expansion coefficients can fully 
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